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ABSTRACT 

We present a study using the Karl G. Jansky Very Large Array (VLA) of 12 CO 
J — 1— emission in three strongly lensed submillimetre-selected galaxies (SMM J16359, 
SMM J14009 and SMM J02399) at z = 2.5-2.9. These galaxies span L m = lO n - 13 L , 
offering an opportunity to compare the interstellar medium of LIRGs and ULIRGs at 
high redshift. We estimate molecular gas masses in the range 2 — 40 x 10 9 M using 
a method that assumes canonical underlying brightness temperature (Tb) ratios for 
star-forming and non-star-forming gas phases and a maximal star-formation efficiency. 
A more simplistic method - using Xco = 0.8 and the measured Tb ratios - yields gas 
masses twice as high. In SMM J14009 we find ^003-2/^001-0 = 0.95±0.12, indicative 
of warm, star- forming gas, possibly influenced by the central active galactic nucleus 
(AGN). We set a gas mass limit of 3a < 6 x 10 8 M for the Lyman-break galaxy, 
A2218 #384, located in the same field as SMM J16359 at z = 2.515. Finally, we use the 
rest-frame -115-GHz free-free flux densities for SMM J14009 and SMM J02399 - mea- 
surements tied directly to the photoionisation rate of massive stars, and made possible 
by VLA's bandwidth - to estimate star- formation rates (SFRs) of 400-600 M yr _1 , 
and to estimate the fraction of Lir due to AGN. 

Key words: galaxies: active — galaxies: high-redshift — galaxies: starburst — sub- 
millimetre — ISM: molecules — galaxies: ISM 
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1 INTRODUCTION 

Observations of the interstellar medium (ISM) in high- 
redshift galaxies can provide insight into the process of star 
formation within these systems. Although the bulk of the 
molecular gas - the fuel for future star formation - is in the 
form of H2 , its lack of a net dipole moment prevents H2 from 
radiating strongly under normal ISM conditions. However, 
this gas can be detected indirectly due to the collisional 
excitation it induces in heteronuclear and near-ubiquitous 
12 CO - the second most abundant molecule in the ISM - 
which emits electromagnetic radiation when its angular mo- 
mentum quantum number changes by %. The low dipole mo- 
ment of 12 CO enables its excitation in regions of low density, 
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^crit ^ 10 2 cm -3 , making it an ideal tracer of the bulk of the 
molecular gas reservoir, allowing us to determine its prop- 
erties in a manner unbiased by the patchy dust extinction 
which can lead to misleading views of velocity fields and 
galaxy morphologies via optical/near-infrared studies. 

Since the earliest detections of 12 CO i n high-redshift 
galax ies (e.g. iBrown &; Vanden Boutl Il99ll : I Solomon et al.l 
1992), much effort has been put into discerning the mass, 
extent and physical properties of the H2 gas in galax- 
ies across cosmic time; this offers a snapshot of their 
evolutionary state and narrows down the possibilities for 
their likely descendents. Until now, available technology 
has limited the study of 12 CO at high redshift to rela- 
tively bright, mid-J lines ( 12 CO J — 3 — 2 and above), 
which are shifted into broad, clean atmospheric windows 
such as that at 3 mm. The excitation requirements of these 
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lines render such observations insensitive to any cool, ex- 
tended gas component, wherein much of the gas is routinely 
found fe-glPapadopoulos fe Seaauistll 19991 : IWeifi et alf 20051 : 
iBothwell et all l2012h . For local galaxies the intrinsically 
fainter J = 1 — transition occupies the 3-mm atmospheric 
window and mid- J lines are more difficult to study, requiring 
high-frequency receivers and excellent atmospheric condi- 
tions. Efforts to conduct comparative studies at low and high 
redshift have therefore been few and far between, and the 
relationship between high-redshift galaxies and their myriad 
possible present-day descendants is not well understood. 

Gravitational amplificatio n - either by chance, line-of- 
sight foreground galaxies (e.g. iRowan-Robinson et al.lll99ll: 



Graham fe Liul Il995l : iNegrello et al.1 l20ld ) or by d eliber- 
ate e xploitation of massive, foreground clusters (e.g. iKneibl 
120101 ) - has provided samples that have formed the basis 
for much of our understanding of the distant starbursts that 
dominate the cosmic far-infrared/submm background. Usu- 
ally buried in the confusion noise, these faint galaxies be- 
come relatively easy to detect and spatially resolve by virtue 
of the magnification they experience. Such an approach led 
to the discovery o f submillimetre-selected galaxies (SMGs - 
ISmail et al.ll 19971 ). to the relatively straightforward identifi- 
cation of their radi o and optical coun terparts, acquisition of 
their redshifts (e.g. Ivison et al. 1998) and their detection in 
mid- J 12 CO (e.g. lFraver et alJ ll998). The r edshift distribu- 
tion of SMGs is broad, centred at z ~ 2.3 ([Chapman et al.l 
l2005h . and their inferred stellar masses (~ 10 10 " 11 M ) 
and star- format ion rates (SFRs; ~ 10 3 M© yr _1 ) make them 
amo ng the most m assive and active galaxies at this epoch 
(e.g. lHainline et aDl2Qllh . 

Here, we follow this familiar pattern. We exploit the 
recent upgrade to the National Radio Astronomy Observa- 
tory's (NRAO's 1 ) Karl G. Jansky Very Large Array (VLA), 
which includes the provision of Ka-band receivers (26.5- 
40 GHz), to study the redshifted 12 CO J = 1-0 emis- 
sion from three lens ed SMGs: ( i ) SMM J 14009+0252 ly - 
ing at z = 2.9344 (|lvison et al.1 120001 : IWeifi "et al.1 l20pdh . 
lensed by Abell 1835 and one of the brightest S MGs detected 
in th e SCUBA Cluster Lens Survey (SCLS; ISmail et al.1 
2002 ) (ii) SMMJ16359+ 6612, a bright SMG a z = 2.5168 
([Kneib et al.ll20o3j2005l ) and lensed into three images by the 
Abell 2218 cluster and (hi) SMM J02399-0136 at z = 2.808 
(jlvison et alJll998h . the first SCUBA galaxy and the bright- 
est of the SCLS galaxies, lensed by Abell 370. We also search 
for 12 CO J=l — emission in (iv) the Lyman-break galaxy 
(LBG) A 2218 #284, located 22.1arcsec from the centre 
of the lensing clus ter in the same field as SMMJ16359 
(jEbbels et alJll996h . 

The optical spectrum of SMMJ02399, coupled with 
its extreme luminosity and emission-line widths (fwzi ~ 
l,000kms _1 ) have long been cited as evidence that this 
source contains a type-2 active galactic nucleus (AGN; 
llvison et a0l998h . SMM J14009 too has long been suspected 
to host a dust-enshrouded AGN on the basis of the strength 
of its radio emission (S ±a ghz — 529±30/xJy) and compact- 
ness (Ivis on et afl koOO. Biggs et al. in prep. ). 

Due to the effects of lensing, we benefit from a boost in 
flux density and apparent size of our targets, surface bright- 
ness being conserved. Excellent lens models exist for Abell 
370, 1835 and 2218 and multi-wavelength imaging has con- 
strained the ages, stellar masses and SFRs of our targets, 



as well as the mass of warm, star- forming gas as traced by, 
e.g., 12 CO J = 3— 2. Observations of another SMG lensed by 
A 1835, SMM J14011+0252 at z = 2.5652, will be discussed 
by Sharon et al. in prep. 

Throughout the paper we use a cosmology with Ho = 
71kms _1 Mpc -1 , Q m = 0.27, Q\ = 0.73. which gives an 
average angular scale of 8.1 kpcarcsec -1 for our sample. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 VLA observations 

Observations of SMMJ14009 and SMMJ16359 were car- 
ried out in blocks of 2-4 hr during excellent weather con- 
ditions between 2010 April and 2011 January (Project IDs: 
AI139 and AI142). Over this Open Shared Risk Observ- 
ing period the available bandwidth from the new Wide- 
band Interferometric Digital ARchitecture (WIDAR) cor- 
relator consisted of two sub-band pairs, each with 64 x 2- 
MHz dual-polarisation channels, and the number of useful 
Ka receivers roughly doubled, from ~ 10 to ~ 20. Due to 
the different redshifts of the targets, different tuning strate- 
gies had to be employed in each case: for SMMJ16359 we 
used both sub-band pairs with the frequency of the red- 
shifted 12 CO J =1 — line (i/ rest = 115.27120256 GHz - 
iMorton fe Noreaul [l994h placed in channel 40 of the lower 
sub-band and a 10-channel overlap between the two sub- 
bands to mitigate the effects of noise in the edge chan- 
nels. This configuration gave velocity coverage and resolu- 
tion of ^2,000 km s _1 and ~18 km s _1 , respectively. For the 
higher-redshift SMM J 14009, the redshifted emission was ac- 
cessible to only one of the two sub-bands and the line was 
placed in its centre providing ^1,300 km s _1 coverage, the 
other sub-band being tuned to z/ Q bs = 32.4 GHz in order to 
constrain the ^130-GHz radio continuum. 

By the time SMM J02399 was observed during 2011 
September through November (Project ID: AT400), the 
OSRO capabilities had been upgraded to include two inde- 
pendently tunable output pairs of eight sub-bands each, with 
128 x 1-MHz dual-polarisation channels per sub-band, giv- 
ing a total available bandwidth of 2,048 MHz. Again, how- 
ever, the redshift of SMM J02399 could be reached by only 
the BD output pair, giving ~20,000kms -1 of coverage and 
^lOkms -1 resolution. The 8 sub-bands of output pair AC 
were tuned continuously at the lowest possible central fre- 
quency, 32.52 GHz, to constrain the radio continuum emis- 
sion from this source. 

For SMMJ16359 we obtained approximately 19 and 
14 hr of useful data in the VLA's D and C configurations; 
for the equatorial target SMMJ 14009 we employed the DnC 
and CnB configurations (with their longer northern arms), 
acquiring 8 and 16 hr, respectively. 10 hr of on-source D 
configuration observations were performed for SMM J02399, 
though an unfortunate technical problem with the correlator 
resulted in only 1/3 of the observed data being recorded. 

Antenna pointing was checked using the C-band 
receivers every 60-90 min, and immediately prior to 
scans of the flux calibrator (3C 286 for SMMJ14009 
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Table 1. VLA sample and observing log. 



Target 


Magnifi- 
-ication 


z 


A 7/„u „ 

(GHz) 


Configu- 
-ration 


Observing dates 


Reference 




SMM J14009 


1.5 


2.9344 


29.23 


DnC 


Sep 25, 29, 30 (2010), 


WeiB et al. 


(^2009) 








-29.36 


CnB 


Jan 21, 24, 26, 27 (2011) 






SMM J16359 


14, 22, 9 a 


2.5168 


32.69 


D 


Apr 12, 14, 20, May 05, 20 (2010), 


Kneib et al. 


(2004) 








-32.93 


C 


Dec 31, Jan 04, 06, 10, 11 (2011) 






SMM J02399 


2.45 


2.808 


29.73 


D 


Sep 27, Oct 26 x2, 


Genzel et al. 


(2003) 








30.69 




Nov 01 x 2 (2011) 







Note: "Images A, B, C, respectively. 



and SMMJ16359; 3C48 for SMMJ02399). We tracked 
amplitude and phase with brief scans of J1642+6856 
(for SMMJ16359), J1354-0206 (for SMMJ14009) and 
J0239-0234 (for SMMJ02399), every 5min. These data 
were also used to determine the bandpass. 

Data were nagged, calibrated and imaged using stan - 
dard AXVS procedures, as outlined bv llvison et al.l (|201lh . 

We constructed data cubes for each source using the 
CLEAN algorithm with 6-MHz (~ 50-kms -1 ) resolution 
by binning the native VLA channels to boost signal-to- 
noise. We used natural weighting and Gaussian uv tapers 
to check for any flux which may have been resolved out in 
the high-resolution images but found no evidence for sig- 
nificant missing flux in any of our sources. The resulting 
synthesised beam sizes are 0.88 x 0.72 arcsec 2 with position 
angles (PAs) 125° for SMM J14009, 1.07x0.72 arcsec 2 at PA 
22° for SMM J16359 and 3.08 x 2.56 arcsec 2 at PA 160° for 
SMM J02399; the integrated maps (Fig. 1) reach down to a 
noise levels of 20, 11 and 20 /iJy beam -1 , respectively. 

Spectra (see Fig. 2) were extracted from the data cubes 
with the AXVS task, ispec, which sums the flux density 
within a box around the source for each plane in the cube. 
Error spectra were produced by extracting five off-source 
spectra and calculating the standard deviations among these 
five sky spectra for each spectral channel. The mean noise 
values per spectral channel for SMMJ14009, SMMJ16359 
and SMM J02399 were 0.09, 0.17 and 0.13 mJy beam" 1 , re- 
spectively. 

Line fluxes are measured by integrating the emission 
spectra across the width of each 12 CO line and then checked 
via the equivalent integrated image. In determining both the 
observational and derived source properties, we treat the 
three images of SMM J16359 as individual sources and then 
co-add the results to derive magnification- weighted averages 
where physically appropriate. 

2.2 Infrared luminosities 

Throughout the paper, we make reference to infrared prop- 
erties of our sources, including calculations of the gas masses 
and SFR derived from rest-frame 8-1,000 /mi luminosities, 
Lir. Here, we present the measurements used to derive these 
quantities. 

W e use archival data from Herschel (Pilb ratt et al.l 
l2010h to measure &50/im , ftso/xm and Ssoo/xm , together with 
published SCUBA flux densities at 750, 850 and 13 5 /im . 
We then exploit the SED library of ICharv fe Elbazl (|200lL 
hereafter CE01) in order to determine Lir for each galaxy. 

For SMMJ02399 we use the mid-IR spectrum of 



iLutz et al.l (l2005h . with data covering A Q bs ^ 750/mi taken 
from llvison et all £l998) . In Fig. [3] we see the excellent agree- 
ment between the data and the best-fit SED template. The 
derived value of L m = (80.7 ± 8. 1) x 10 11 L Q is 30 per cent 
lower than the value reported in llvison et all fl998). 

For SMM J16359 we use publish ed flux den s ities f or the 
three images at 24 and 70 fim from lKneib et al.1 (2004). The 
three images are resolved by Herschel up to 350 /im but be- 
come blended at longer wavelengths. Where the three im- 
ages are blended, we divide up the observed flux in pro- 
portion to the known magnification factors. We derive an 
estimate of Lir for each of the three images of this source 
and combine our three estimates to derive the magnification- 
weighted mean Lir = (6.4 ± 1.0) x 10 11 L Q . This va lue is 
consistent with the results of Finkel stein et al. I (|201lh . who 
determined Lir = 7.0 x 10 11 L (for image B) using a similar 
SED-fitting method. 

Similarly, fo r SMM J 14009 we rel y on the 24-/im mea- 
surement from iHempel et al. (2008) to supplement the 
archival Herschel and SCUBA data, and compute a lensing- 
corrected L m = (160 ± 24) x 10 11 L . 

The best-fit CE01 SEDs for each galaxy are shown with 
the corresponding photometry in Fig. 3. 



2.3 Radio continuum data 

The VLA tuning strategies employed to observe 12 CO J — 
1-0 in for SMMJ14009 and SMMJ02399 enabled us to 
obtain detections of the ~ 115-GHz continuum emission 
from these galaxies. These data were reduced in parallel 
with and according to the same recipe as the 12 CO spectral- 
line data. We do not have continuum data at this frequency 
for SMM J16359, owing to the different tuning strategy em- 
ployed at the VLA for this source. We detect continuum 
flux densities of S 3 2GHz = 57 ± 25/xJy for SMM J02399 
and &ghz = 42 ± 20/iJy for SMM J14009. 1.4-GHz flux 
densities of Si a gh z — 5 2 6 =b 5 and 529 =b 30 ^ J y were de- 
termined in llvison et all (|l998h and llvison et aD (|2000h for 
SMM J02399 and SMM J14009, respectively. The VLA sci- 
ence archive was queried for continuum data at 5- and 8- GHz 
as well, with flux densities reported in Table. 4. The radio 
continuum data are plotted with the IR SEDs in Fig. 3. 
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Figure 1. Integrated 12 CO J = 1-0 maps of our three targets, SMMJ16359 (10 arcsec/side), SMMJ14009 (10 arcsec/side) and 
S MM J02399 (15 arcs ec/side). For SMM J16359 - a triply imaged source - we show only the brightest image, SMM J16359-B as denned 
bv lKneib et al 1 d2004h . Main images: 12 CO J = 1—0 contours spaced at —3, 3, 4, 5...xcr, set on top of colour representations of the same 
data. The VLA beam used to produce each map is shown in the bottom-right corner. Insets: 12 CO J = 1 — contours spaced at —3, 3, 
A/2X3...XCT and in steps of y/2 x a thereafter, set on top of three-colour HST images of the same field of view, in order to explore the 
relationship between the SMG and any optical counterpart there may be. For SMM J14009 we detect no optical counterpart to the 12 CO 
emission. 




Figure 2. 12 CO J = 1 — spectra (yellow histograms) and reference J = 3 — 2 spectra (dotted lines) for our three targets. The left-hand 
scale represents J = 1 — flux density, and the right-hand side shows the corresponding J = 3-2 flux density, scaled by 9 _1 so as to 
be on the same Rayleigh- Jeans scale. The zero-LSR redshifts used to determine the relative velocities for each source are quoted 
in the figure, and typical error bars are shown in the top-left-hand corner of each plot. The spectrum shown for SMMJ16359 is the 
magn ification- weighted mean of the spectra of the three images of this source, SMMJ16359-A, -B and -C as defined in iKneib et al.l 



3 ANALYSIS AND RESULTS 

3.1 12 CO morphologies and spectra 

We summarise the source properties in Table 2. All sources 
are detected in 12 CO J —1 — with signal-to- noise of 13 
for SMMJ16359-B (the brightest of three images of this 
multiply- lensed source), 17 for SMMJ14009 and 7.5 for 
SMMJ02399. The flux- weighted 12 CO J =1-0 redshifts - 
summarised in Table 2 - are derived by fitting a single Gaus- 
sian component to each 12 CO J = 1 — spectrum and are 
consistent with the reference 12 CO J = 3— 2 redshifts to four 
significant figures. The spectrum shown for SMM J16359 is 
a magnification-weighted average of the spectra from each 
image of this source. The uncertainties in the spectral fits 
translate to velocity uncertainties that are typically of order 
AV = 10-15 km s" 1 , except in the case of SMMJ02399 
where the comparatively low signal-to-noise hampers our 
ability to constrain zco and results in AV ~ 50kms _1 . 
For each source we see broad lines with full width at zero 
intensity (fwzi) of ~500-l,000 km s _1 . 

We compare the line profiles and maps between J = 1—0 
and J = 3— 2 observations. The astrometry is consistent. For 



SMMJ16359 and SMMJ14009 we see good agreement be- 
tween the low- J and mid- J 12 CO line profiles. In the former, 
we detect double-peaked J= 1—0 line profiles fo r each of the 
three images, separated by ^300 km s -1 , as did lKneib et al.l 
(|2005h in their study of the J = 3- 2 emission. SMM J02399 
displays a significant offset in velocity between the measured 
J — 1—0 emission and the J = 3— 2 spectrum of iGenzel et al.l 
(2003), ^100 km s , roughly twice that which could be as- 
cribed to uncertainties in the Gaussian fits. However, the 
J = 1 — and J = 3 — 2-derived redshifts are consistent to 
four significant figures (z = 2.808). It is difficult to envis- 
age any plausible physical origin for these offsets: in an ef- 
fort to track down any instrumental contribution to the fre- 
quency offsets, continuous- wave tones were injected into the 
Ka-band receivers at frequencies close to those used for the 
science observations, but no offsets were observed between 
the input and observed tones (Emmanuel Momjian, private 
communication) . 

We find a weighted mean velocity width ratio of 
<7(i_o)/(3_2) = 1.3±0.1 for th e three images of SM M J16359, 
consistent with the findings of llvison et al.1 (|201lh for a sam- 
ple of four unlensed SMGs. For the two more AGN-like 
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Figure 3. Observed-frame infrared SEDs for all of our sources, uncorrected for lensing, including data {symbols: see legend) and the 
corresponding best-fit CE01 SED. We include only data up to A Q b s = 1350 in the fit. For clarity, we plot the SED of only one image 
of SMM J16359, and plot the best-fit CE01 template out to radio wavelengths, demonstrating the "canonical starburst" radio SED. For 
the other two galaxies, we truncate the CE01 SED at 2,000 /mi and plot our modelled radio SED, constrained by our measurement of the 
Ka-band continuum flux and previously published results at longer wavelengths ( 32.31 33. 5|) and interpret the steepness of the resulting 
spectra as being consistent with the known AGN hosted within these galaxies. 



Table 2. Observed properties. 



Target ^CO(l-O) ^CO(l-O) ^CO(3-2) r 3-2/l-0 °CO(1-0) cr CO(l-0)/CO(3-2) FWZI R 

(Jykms -1 ) (Jykms -1 ) (kms -1 ) (kms -1 ) (kms -1 ) kpc 



SMM J14009 


2.9332 ± 0.0001 


0.31 ±0.02 


2.70 ±0.30 


0.95 ±0.12 


175 ± 10 


1.06 ±0.08 


850 


1.2 ±0.4 


SMM J02399 


2.8083 ± 0.0007 


0.60 ±0.12 


3.10 ±0.40 


0.58 ±0.14 


350 ± 90 


1.05 ±0.26 


1300 


26 ±4 


SMM J16359-A 


2.5174 ±0.0002 


0.22 ±0.04 


1.67 ±0.13 


0.83 ±0.15 


195 ± 15 


1.10 ±0.09 


710 


< 1 


SMM J16359-B 


2.5173 ±0.0001 


0.40 ±0.04 


2.50 ±0.12 


0.70 ±0.22 


190 ± 10 


1.26 ±0.07 


780 


0.6 ±0.3 


SMM J16359-C 


2.5179 ±0.0002 


0.30 ±0.09 


1.58 ±0.17 


0.59 ±0.20 


305 ± 15 


1.81 ±0.09 


670 


4±2 



galaxies, SMM J14009 and SMM J02399, we detect velocity 
width ratios consistent with unity (<7(i_o)/(3-2) = 1.06±0.08 
and 1.05 ± 0.26 respectively). 



3.2 Line luminosities and derived total gas masses 

Traditionally, attempts to determine the molecular gas mass 
of distant galaxies via 12 CO emission have been subject to 
two main sources of uncertainty: (1) where mid- J lines are 
observed, a brightness temperature ratio, r( 3 _ 2 )/(i-o) = 1 
has typically been assumed to convert the measured 12 CO 
J = 3 — 2 flux to an inferred J = 1 — flux, and (2) the 



"X- f actor" , Xco = M(H 2 )/L / CO(1 _ 0) (jPownes & Solomon! 
1998j) which is used to convert the observed or inferred 
12 CO J — 1 — luminosity to a total molecul ar gas mass. 
As a number of recent papers have shown (e.g. lHarris et al.l 
l2010l : llvison et all l201ll : iDanielson et al.1 boill . etc.), ob- 
served r( 3 _2)/(i-o) ratios in SMGs are typically well below 
unity. We can thus improve our gas mass estimates by ob- 
serving the 12 CO line which traces the bulk of the molecular 
gas reservoir, but we must still worry about the order of mag- 
nitude range spanned by Xco depending on the assumed 
state of the underlying ISM: Xco ^ 5 in optically-thick gi- 
ant molecular clouds (GMCs) where the gas is reducible to 
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Table 3. Lensing-corrected properties 



Target 


^CO(l-O) 


-^gas , ol 


m CQ/SF 


^IR 


A^gas,Ill 


M dyn 




(xlO 9 Kkms-ipc 2 ) 


(xlO 9 M ) 




(xlO 11 L ) 


(xlO 9 M ) 


(xlO 9 M ) 


SMM J14009 


83 ±6 


67 ±4 


0.09 ± 0.02 


160 ± 16 


35 ±7 


18 ±4 


SMM J02399 


96.6 ±1.9 


77 ± 15 


1.5 ±0.5 


81 ±8 


41 ± 14 


1560 ± 550 


SMM J16359-A 


4.9 ±0.8 


3.9 ±0.6 


0.33 ±0.09 


7.1 ±0.7 


1.3 ±0.4 


< 19 


SMM J16359-B 


5.5 ±0.6 


4.4 ±0.5 


0.8 ±0.3 


6.2 ±0.6 


2.8 ±1.3 


11 ±4 


SMM J16359-C 


10 ±3 


8±2 


1.4 ±0.7 


5.7 ±0.6 


1.5 ±0.7 


180 ± 70 


A2218 #384 


3cr < 0.8 


3cr < 0.6 











Note: a Calculations for Lir discussed in £12.21 All luminosities have been corrected for the known magnification factors. 



ensembles of self- gravitating units, whereas Xco ~ 0.8 in 
the more extreme environments of local ultra luminous in- 
frared galaxies (ULIRGs), where the ISM forms a smooth, 
continuous medium. This latter value is typically adopted 
for SMGs. 

The intrinsic 12 CO J — 1 — line luminosity of an un- 
lensed galaxy is the observed line flux integrated across 
the surface area and velocity range of the line: L co = 



Iav I a 



S u dA dV K km s 1 pc 2 , or: 



^co(i-o) = 3.25 x 10 7 



D 2 (Mpc) 



AV 



1 + z 

S u dis 



GHz 



(l) 



Jy km s 



pc 



Having made an appropriate correction for magnifi- 
cation effects, we then invok e the commonly-used CO- H2 
gas mass conversion factor of iDownes &; Solomonl (|l998l ) of 
Xco ~ 0.8 M (Kkm s _1 pc 2 ) -1 to derive total gas masses 
(Table 3). 

Combin ing our 12 C O J — 1 —0 data with existi ng J — 3—2 
data from iKneib et al.1 ([20051, SMMJ1635 9), IWeifi et al.1 
(|2009l . SMMJ14009) and lGenzel eTal] (|200l SMMJ02399) 
enables us to derive a mean magnification-weighted Tb ra- 
tio of r( 3 _2)/(i-o) = 0.72 ± 0.12 for the three images of 
SMM J16359, r (3 _2)/(i-o) = 0.95±0.25 for SMM J14009 and 
r(3-2)/(i-o) = 0-58 ± 0.14 for SMM J02399: SMGs on aver- 
age experience sub-thermal excitation. We can then compare 
the gas mass derived above (M g as ,a) with a gas ma ss derived 
according to the prescription of Iviso n et al 1 (1201 ih . denoted 
Mgas.in, in which a lower limit for the mass of actively 
star- forming gas (Msf) is computed on the assumption that 
SMGs form stars at an appreciable fraction of Eddington- 
Limited star- format ion efficiency (SFE max = Lir/Msf ^ 
500 L M^ 1 where Lir i s the lensing-corrected infrared lu- 
minosity: IScovillell2004h . For such maximally- efficient star 
formation, a measurement of Lir therefore corresponds to a 
determination of Msf- Using canonical brightness tempera- 
ture ratios for cold, quiescent (^[^?2)/(i-o) = 0.3) an d warm, 

star- forming (^[f^2)/(i-o) = & as ' their relative fractions 
in a sub-thermally-excited reservoir can be computed as 



^CQ/SF — 



(SF) 

r (3-2)/(l-0) ~ r (3-2)/(l-0) 
,(CQ) 



(2) 



^(3-2)/(l-0) - > ( 3 -2)/(l-0) 

and the total molecular gas mass is calculated as 



M gas = M SF (1 + mcQ/sr) (3) 
Computing M gas in this way provides an indepen- 



dent check on the total gas mass, based on different as- 
sumptions from those that underpin the X-factor method. 
The magnification- weighted mean value of M gas ,iii for 
SMM J16359 is (2.1 ±0.6) x 10 9 M , and for SMM J02399 we 
detect (41.0 ± 14.4) x 10 9 M of gas using this method. For 
SMM J14009 we calculate M gas ,iii = (34.9 ± 7.3) x 10 9 M . 
All three values of M ga s,iii are ~ 2x lower than the equiv- 
alent values calculated via Xco, which would imply Xco ~ 
0.4. We remind the reader, however, that M ga s,ni is calcu- 
lated on the assumption that SMGs have maximal SFE: if 
the SFE is less than 500 L Mq 1 , then a higher mass of star- 
forming gas is required for a given Lir, raising the total gas 
mass in the system accordingly. 

We can then compare these gas masses to the dynamical 
masses, assuming Md yn = 2.1 Ra co ^ 1 _ ^ ) /G dTacconi et al.l 
2008). Radii for our sources are derived by fitting a two- 
dimensional Gaussian to each integrated 12 CO J — 1—0 map. 
This estimate assumes that the gas emission traces a vit- 
alised potential well. Using the observed values of oco(i-o) 
reported in Table 2, we determine dynamical masses of 
(18±4) x 10 9 M , (11±4) x 10 9 M and (1560±550) x 10 9 M 
for SMMJ14009, SMMJ16359-B and SMMJ02399, imply- 
ing firm upper-limits of Xco < 0.3, Xco < 2.8 and Xco < 
24, respectively. The dynamical masses derived from the two 
brighter images of SMMJ16359 are ^3x higher than the 
best-estimate of the gas mass derived via X co = 0.8, consis- 
tent with the sample of Ivis on et al. I (|201lh . The extremely 
compact radius quoted for SMM J 14009 is a result of the 
narrowness of the 2D Gaussian that best fits the integrated 
map, itself a result of the source being only marginally re- 
solved on North- South baselines and completely unresolved 
on East- West baselines. The gas in SMMJ02399 was re- 
ported as h a ving a lensing-corrected extent of ^25kpc in 
llvison et al.l (|2010h in their C-configuration observations, 
corroborating the equally large extent of 26 ± 4 kpc derived 
from our D-configuration observations, as reported in Ta- 
ble 2. The colossal value of Md yn therefore likely implies that 
the assumption the gas traces a virialised potential well is 
false. 



3.3 Gas excitation: observed brightness 
temperature ratios 

Having worked with canonical brightness temperature ratios 
to estimate M gas , free from any assumptions about Xco, we 
now compare our 12 CO J — 1 — spectra and the literature 
J — 3 — 2 spectra in order to determine the true excitation 
conditions of the gas in our samples of galaxies. 
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The steepness of the observed radio continuum spec- 
trum of SMM J02399 f ffPl) is consistent with the claim that 
this source harbours a powerful AGN. Such an AGN might 
be expected to heat the surround ing gas via bombardment 
with X-rays (|Meiierink et alj20Qq ) or through shock-heating 
as the po werful radio jet drives into the surrounding m olec- 
ular gas ( Pa padopoulos et al.ll201ll : llvison et al.l 12012). yet 
the Tb ratio detected for this source r( 3 _ 2 )/(i-o) = 0.58 ± 
0.1 4 is entirely consist ent with that of less AGN- like SMGs. 
As llvison et all (|2010h observe, however, SMM J02399 ap- 
pears to consist of three or more interacting galaxies em- 
bedded within a larger reservoir of cold molecular gas. In 
particular, the massive obscured starburst which gives rise 
to the extreme Lir of this source is spatially offset from the 
broad-absorption-line quasar. 

Our improved sampling of the radio SED of 
SMMJ14009 allows us to determine that this source also 
has a steep radio continuum spectrum ( fc|2. 3p , corroborat- 
ing its classification as an AGN. Additionally, SMM J14009 
demonstrates a high Tb ratio, r( 3 _ 2 )/(i-o) = 0.95 =b 0.12, 
which indicates the presence of a warm gas phase given the 
excitation requirements of the 12 CO J = 3 — 2 transition, 
n crit - 10 4 cm- 3 , E u /k B ~ 33k. 

For SMM J16359, the magnification- weighted mean Tb 
ratio, r( 3 _2)/(i-o) = 0.72 ± 0.12, is entirely consistent with 
that of the SMG population ([Harris et alj|201ol : llvison et al.l 

[mil). 

In addition to calculating source- averaged Tb ratios, our 
new data give us the spatial and spectral resolution to assess 
r( 3 _2)/(i-o) as functions of position and velocity. In Fig.[3]we 
illustrate the variation in gas excitation across the observed 
line profiles for the three sources. The results are only con- 
sidered to be robust within the central ±400 km s -1 of each 
of the line profiles due to deteriorating signal/noise at the 
edges of the line profiles; nevertheless, in SMMJ16359 we 
see clear structure in the plot of r( 3 _ 2 )/(i-o) versus velocity, 
with r( 3 _2)/(i-o) ~ 1 for regions in the centre of the velocity- 
space, rolling off to sub- thermal values for the high- velocity 
wings of the line. These variations in the velocity fields of 
each galaxy may point to spatial variations in r( 3 _ 2 )/(i-o), 
and are consistent with the observed values of 0"(3_2)/(i-o) • 

We attempt investigate the velocity structure of our 
galaxies by constructing position- velocity (PV) diagrams for 
each source, extracting a slice through the (right-ascension, 
declination, velocity) data cube along the major axis of the 
galaxy. The resultant PV diagrams are shown in Fig.[5j how- 
ever their low signal-to-noise prohibits us from modelling the 
emission. 

3.4 Infrared- derived star- for mat ion rates 

We present our infrared SEDs in Fig. [3] and use the in- 
frare d luminosities repor ted in Table 3 with the t runcation 
from iGenzel et aD (|2003h of the iKennicuttl (1 19981, here after 
K98) conversion factor for a 1-100 M ISalpeterl ((,1955) ini- 
tial mass function (IMF) to compute SFRir: 

SFRir _ Lir . . 

1 Moyr- 1 1.2 x 10 10 L [ } 

In the case of SMMJ02399, Lir is thought to arise 
roug hly equally from the AGN and starburst components 
(e.g. IValiante et al.1 l2007h . Adopting this AGN fraction 
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Figure 4. Top - bottom: Plots of Tb ratio as a function of 
velocity for the composite SMMJ16359 spectrum, SMM J14009 
and SMMJ02399. Due to poor signal/noise at the extremities 
of the line profiles, we only consider data within the central 
±400 kms -1 to be robust enough to draw any conclusions from. 
The composite spectrum for SMM J16359 is computed by taking 
the magnification-weighted mean of the individual spectra from 
images A, B and C of this triply imaged source. The observed 
velocity offset between 12 CO J = 1 — and J = 3 — 2 spectra in 
SMM J02399 manifests itself as a super-thermal Tb ratios in the 
envelope of the velocity profile. Typical uncertainties are denoted 
by the error bar in the top-left of each plot. 

for both SMMJ02399 and SMMJ14009, Equation H yields 
SFRir = 340 ±40 and 670 ±70 M yr" 1 , respectively, where 
the uncertainties are purely statistical (due to Lir), and thus 
do not include the uncertainty in the K98 conversion factor. 

3.5 An independent probe of the SFR via 
optically-thin free-free emission 

The SED longward of the FIR dust bump in galaxies is char- 
acterised by three main components: (i) the Ray leigh- Jeans 
tail of the infrared/sub-mm dust spectrum, (ii) thermalised 
free- free emission arising from ionised Hn regions, which 
traces the massive/young stars (> 5 M©) capable of pho- 
toionising, and (iii) synchrotron radiation ar ising from the 
electrons accelerated by supernova remnants (iCondoni ri992. 
hereafter C92). Sometimes (iv) synchrotron radiation aris- 
ing from plasma accelerated by an active AGN may also 
contribute. 

Our data and the new estimates of Lir provide us with 
the means to determine the SFR via two independent meth- 
ods: (i) via Lir between 1-100 M©yr _1 (K98) and (ii) via 
the optically-thin free-free emission, probing SFR(^ 5M©) 
(C92). 

We determine SFR ra di via Equation [5] where Lff is the 
free- free luminosity density at frequency z/, defined such that 
L = And 2 S& for lensing-corrected flux density, Ss. The free- 
free luminosity is derived by collating the radio continuum 
data summarised in Table 4, and fitting the relative contri- 
butions to the observed SED arising from dust (an extrapo- 
lation of the Ray leigh- Jeans side of the best-fitting CE01 in- 
frared template), thermalised free-free emission (Ss oc z/ -0 " 1 : 
C92) and synchrotron emission (S'sync oc is~ async ). The free 
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Figure 5. Clockwise from top left: Position- velocity diagrams for SMM J16359-A, SMM J16359-B, SMM J16359-C, SMMJ02399 and 
SMM J14009. 



parameters are the flux-scaling of the synchrotron and free- 
free emission, and also the spectral index, a sy nc, of the syn- 
chrotron component. 



L ff 



WHz- 



5.5xl(r 



GHz 



SFR(M ^5M G 
Moyr- 1 



(5) 



Using the 1.4-, 5- and 8-GHz data available for our two 
galaxies, we measure spectral indices, a (where S u oc 
of -0.96 ± 0.06 and -0.98 ± 0.14 for SMMJ14009 and 
SMM J02399, respectively, a little s teeper than the m ean 
value of -0.75 ± 0.06 measured by llbar et al.1 (|2010h be- 
tween 610 MHz and 1.4 GHz for 44 radio-detected SMGs in 
the Lockman Hole, and steeper too than the CE01 template 
that best fits each galaxy's mid-IR-to-submm SED. 

The contributions to S32 GHz arising due to dust, syn- 
chrotron and free- free emission are presented in Table 5. 
We determine SFR(M ^ 5 M ) = 290 ± 85 Moyr -1 for 
SMMJ14009 and 200 ± 20 Moyr -1 for SMM J02399, and 
extrapolate down t o 1 Mq. Accou nting for the unseen low- 
mass stars using a ISalpeterl (|1955l ) IMF, we find 1-100 M© 
SFR radio of 630 ± 190 Moyr" 1 and 430 ± 50 Moyr" 1 for 
SMMJ14009 and SMMJ02399, respectively. In common 
with the preceding section, these uncertainties are statisti- 
cal, and do not reflect the uncertainty in the C92 conversion 
factor. 



3.6 The Lyman-break galaxy, Abell 2218 #384 

Returning to the observations for SMM J16359, we searched 
for 12 CO J = 1-0 in the LBG at (R.A. 16:35:49.4, Dec. 



66:13:07 J2000: lEbbels et al.|[l996h located behind the same 
cluster. We do not detec t this source, but set an upper limit 
using S ^ 3cr n V AvAw (|Seaquist et al . 1995), where Av is 
the velocity resolution of t he spectrum, Aw ~ 250 km s" 1 
is a typical LBG line width (|Greve &; Sommer-Lar sen 2008) 
and the noise in the primary beam-corrected VLA spectrum 
is a n . We de-boost this by ~ 16 x to accoun t for gravitational 
lensing by the cluster (|Ebbels et alJll996h and thus derive 
upper limits on the 12 CO J —1 — luminosity of < 0.8 x 
10 9 Jy kms -1 pc 2 and on the gas mass of < 6 x 10 8 M©, 
respectively, where we have assumed Xqo = 0.8, supposedly 
appropriate for star- forming galaxies. 

This gas mass limit is comparable to gas masses de- 
tecte d in two similarly high-redshift LBGs by Rie chers et al.l 



2010h . who determined M gas = (9.3 ± 1.6) x 10 8 M and 
M g as = (4.6±1.1) x 10 8 M for the same A C o, in the Cosmic 
Eye (at z = 3.074) and MS 1512-cB58 (at z = 2.727). 



4 DISCUSSION AND CONCLUSIONS 

4.1 The physical conditions of the molecular gas 

We have presented high-resolution VLA imaging of 12 CO 
J — 1 — in three bright, lensed SMGs. Two of our sources 
harbour known AGN - SMMJ14009 and SMM J02399 - 
an d are brighte r (Lir ~ 10 13 L©) than the unlensed sample 
of llvison et al.l (|201lh . while the other source - the triply 
imaged SMMJ16359 - is an order of magnitude fainter 
(Lir rsj 6 x 10 11 Lq, observed as being equally bright due to 
the high magnification afforded by its lensing cluster, Abell 



© 2012 RAS, MNRAS OOO.fTUTuI 



VLA imaging of ionised and molecular gas in SMGs 9 



Table 4. Observed radio properties 



Target 


Si. 4 GHz 


#5 GHz 


#8 GHz 


#32 GHz 




name 


(MJy) 


(A*Jy) 


(A*Jy) 


(^Jy) 




SMM J14009 


529 a ± 30 


151 ± 17 


91 ± 18 


42 ± 20 


0.96 ± 0.06 


SMM J02399 


526 b ± 50 


130 ± 22 


N/A 


57 ±25 


0.98 ± 0.14 



Notes: " llvison et all d2000h : ^Ivison et ahl dl998h . 



Table 5. Radio source derived properties 



Target 
name 


^Ka,sync 

(MJy) 


C^sync 


^Ka^ust 

(A*Jy) 


Ska,ff 

(A*Jy) 


SFR a 

radio 

Moyr- 1 


SFR^ 
Moyr" 1 


SMM J14009 
SMM J02399 


22 ± 2 
16 ± 4 


1.01 db 0.09 
1.09 ± 0.13 


8 ± 1 
18 ± 1 


12 ± 5 
24 ± 6 


630 ± 190 
430 ± 50 


670 ± 70 
340 ± 40 



Notes: a SFRs calculated for 1-100 Mq Salpeter IMF; b SFRiR calculated on the assumption that the AGN and star formation 

contribute equally to Lir. 



2218), and does not appear to be dominated by AGN emis- 
sion. 

Comparing integrated maps of 12 CO J= 1-0 and J = 3- 
2 emission, we compute mean brightness temperature ratios 
of less than unity and velocity width ratios a( ^_ 2 )/(i-o) Z 
1, bot h results consistent with the findings of llvison et al.l 
(2011): SMGs on average experience sub-thermal excitation 
and exhibit broad 12 CO J=l — line widths. 

We derive Lir for our galaxies by collating litera- 
ture data, and then choosing from the SED library of 
ICharv fe Elbazl (|200lh the best-fitting SED, which we in- 
tegrate between rest-frame 8-1,000 /xm. Next, we compute 
molecular gas masses via two methods: first by converting 
Leo to Mgas via an X- factor of Xco = M(H2)/I/ C o(i-o) 
and then, as an independent check, using an approach that 
assumes SMGs form stars with Eddington-limited SFE, 
which gives us a lower-limit to the molecular gas mass. For 
each of our sources, the latter method yields a gas mass 
~ 2x lower than that derived via Aco- 



4.2 Independent estimates of the SFR, and the 
role played by AGN 

We compare SFRs derived via two independent diagnostics: 
(i) SFRir, as computed via Lir (K98) and (ii) SFR ra dio, 
computed from Lff, the optically-thin free- free emission 
(C92) (Table 5). Both of these diagnostics are unaffected by 
interstellar extinction, however the former method is liable 
to over-estimate the SFR if a signific ant fraction of Lir a rises 
due to heating of dust by an AGN (|Murphv et al.ll201lh . In 
contrast, the decomposition of the radio SED into dust, syn- 
chrotron and free- free contributions provides us with an SFR 
diagnostic that is directly related to the photoionisation rate 
of young, massive stars, relatively unaffected by AGN. 

To derive SFRir in §3.4, we assumed that ~50 per cent 
of Lir was due to star form ation, with the rest due to AGN 
heating of the dust (e.g . iFraver et al. | Il998l : iGreve et al.1 
2005; IValiante et all 1200/T ). Here, we use our new measure- 
ments of free-free emission - which is unaffected by the pres- 
ence of an AGN - to test this assumption. Using Equation 
H without an AGN correction, SMMJ02399's SFRiR, upP er 
is 670 ± 70 M yr -1 . This is 35 ± 10 per cent higher than 



SFR ra dio for a 1-100- M© Salpeter IMF, implying a -35 
per cent contribution to Lir from its AGN. Similarly for 
SMM J14009, SFRiR, U p P er = 1,340 ± 140 Moyr" 1 , indicat- 
ing that some 55 =b 15 per cent of Lir arises due to the AGN. 

Recognising that the spatial extent traced by 12 CO 
J — 1 — for each of our galaxies is almost certainly larger 
than the region in which star formation is actively taking 
place, we can use our new free- free- derived SFRs to place 
lower-limits on the global star- format ion rate surface den- 
sity, Ssfr = 0.7 ± 0.1 M yr- 1 kpc~ 2 for SMM J02399 and 
Ssfr = 450 ± 200 M vr -1 kpc~ 2 for SMM J14009, the for- 
mer value again attributable to the colossal extent of the 
12 CO J=l-0 reservoir in SMMJ02399. We note that the 
peak Ssfr will be significantly higher than these galaxy- 
averaged values. 

At present, our measurement of SFR ra dio - and hence 
our estimates of the AGN fraction and Ssfr - are only 
as good as the extrapolation of the best-fit SED to the 
IR/submm data, via which we determine the dust contribu- 
tion to the z/ bs °° 11 5- GHz flux density, itself a relatively low 
signal-to- noise measurement. However, these uncertainties 
can be reduced significantly via better sampling of the SED 
beyond the wavelengths probed by Herschel. When com- 
plete, the VLA will offer 8 GHz of simultaneous bandwidth 
across 64 independently tunable sub-band pairs. Together 
with the Atacama Large Millimeter Array, it will be pos- 
sible to measure the 350 /xm-20 cm SEDs of these galaxies 
with an order of magnitude better sensitivity than currently 
possible. 
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